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Sunspot umbra atmosphere from full Stokes inversion 
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Abstract. Sunspots are prominent manifestations of the solar cycle and 
provide key constraints for understanding its operation. Also, knowing internal 
structure of sunspots allows us to gain insights on the energy transport in strong 
magnetic fields and, thus, on the processes inside the convection zone, where 
solar magnetic fields are generated and amplified before emerging at the surface 
on various scales, even during solar minima. In this paper, we present results of a 
spectropolarimetric analysis of a sunspot observed during the declining phase of 
the solar cycle 23. By inversion of full Stokes spectra observed in several spectral 
regions in the optical at the THEMIS facility we infer the height dependence 
of physical quantities such as the temperature and the magnetic field strength 
for different sunspot regions. The simultaneous use of atomic (Fe I 5250.2 and 
5250.6 A) and highly temperature sensitive molecular (TiO 7055 A and MgH 
5200 A) lines allow us to improve a model of the sunspot umbra. 



1 Introduction 

Inversion of spectropolarimetric data becomes a standard tool for re covering; the 



vertical structure of solar magnetic features, including sunspots (e. g. jMathew et al 



20Ql : lBeckl[2008l ). However, in most cases the inversion is limited to atomic lines 



formed in a narrow range of heights in the solar atmosphere. Simple diatomic 
molecules like MgH and TiO are sensitive to both the magnetic field and the 
temperature and provide inform ation on different atmospheri c layers of the solar 



)spnenc laj 

phot osphere than atomic lines (Be rdyugina fc Solankill2002l : iBerdvurina et al. 



120031 ). They can be observed in the sunspot umbra at typical temperatures be 
low 4500 K. Combining spectropolarimetry of atomic and molecular lines thus 
provides a very sensitive tool for detecting and probing magnetic structures 
i n the sunspot umbr a, which so far was done for infrared Fe I a nd OH lines 
dMathew et all l2003h and for optical Fe I, MgH, and TiO lines dAfram et al ' 



2006). In this paper, we further employ the advantage of this approach. 
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2 Observations 



The sunspot NOAA 10667 was observed at the THEMIS facility ( Tenerife) din- 
ing se veral consecutive days in September 2004 ( Afram et al.ll2006l : lArnaud et al. 



2006). Full Stokes profiles were recorded simultaneously in five spectral regions. 
From the available data three spectral windows were chosen for the present anal- 
ysis: Fe I lines at 5250.2 A and 5250.6 A, TiO at 7055 A, and MgH at 5200 A. 
As indicated by contribution functions shown in Fig. [3J the chosen atomic and 
molecular lines carry information on different layers in the sunspot atmosphere, 
where they constrain our model. In addition, their different Lande factors allow 
us to deter mine both th e strength of the magnetic field and the magnetic fill- 
ing factor (jStenflo Ill973l ). This sunspot was also imaged during a coordinated 
program at the SST (La Pal ma) with a narrow-band filte r centered at the TiO 
7055 A band, as designed by lBerger fe Berdyuginal ( 20031 ) (see Fig. [1]). 




Figure 1. Left: Image of the sunspot NOAA 10667 obtained in the TiO 
7055 A band at the SST. The resolution and position of the THEMIS data 
are indicated by the grid. The top and right hand side scales give the angular 
size of the sunspot which is about 15x25 arcsec 2 . For each pixel the full 
Stokes vector was recorded. Right: The results of the inversion. These maps 
show the temperature and magnetic field stratifications at log T5000 = — 1 for 
the cool second magnetic component. 



3 Inversions 

Stokes parameters I/I c , Q/I c , U/I c and V/I c in the selected spectral regions were 
inverted for each pixel indepe ndently employing the code SPINOR fjFrutiger et al 
20001 ; iBerdyugina et al.l l2003). This code assumes LTE conditions and solves the 
Unno-Rachkovsky radiative transfer equations using response functions. Start- 
ing from an initial guess the Stokes profiles are fitted iteratively to the obser- 
vational data by varying parameters of the model atmosphere. In Fig. E] an 
example is shown when temperature and magnetic field strength are iterated as 
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height dependent free parameters at five reference points log T5000 — — 4, —3, 
-2, -1, and 0. 

After extensive tests with different atmospheric models, a three component 
atmosphere was found to give the best fit results. The first component is a 
standard quiet sun atmosphere accounting for straylight, whereas the second 
and third component represent magnetic parts of the sunspot atmosphere. All 
spectral windows and thus all atomic and molecular lines were inverted simul- 
taneously. 



1.04 
0.85 
0.66 
0.47 
0.28 
8.30 


-8.30 









V 






/K/v 







-0.40 0.40 
5250.5 


X-X [A] 
7054.4 


7055.1 


-0.40 0.40 
5199.6 





-0.40 0.40 X-X [A] 

5250.5 7054.4 7055.1 



-0.40 0.40 
5199.6 



-0.40 0.40 X-X [A] 




-0.40 0.40 l-X [A] 



Figure 2. The best fit to a typical set of Stokes profiles for one umbral pixel 
located at (7, 53) according to the grid shown in Fig. [TJ The dotted line is 
observational data, the solid line is the synthetic spectrum found by inversion. 



4 Results 



Even though the inversion was applied to the whole set of data, here we discuss 
results for the sunspot umbra only. As indicated above, the resulting atmo- 
sphere leading to best fits consists of three components. The overall magnetic 
filling factor is about 0.92 splitted almost equally between the two magnetic 
components. These components have different temperature gradients near the 
bottom of the atmosphere (see Fig. [3D and most probably indicate the presence 
of unresolved umbral structures such as umbral dots imbedded into a diffuse 
dark umbra. 

We compare our results with models of the dark core umbra by lMaltby et al.1 
f|!986h . Those are also shown in Fig. [3] and correspond to early (E) and late (L) 
phases of the solar cycle or may be interpreted as having the properties of dark 
and bright umbra, respectively. Clearly, our umbral models are warmer than 
those of Maltby. This is in agreement with a systematic inc rease of the umbral 
tempe rature during at least the last decade as reported by iPenn fc Livingston 
f|2006h . 
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Figure 3. A comparison of different umbra models. Average umbral tem- 
perature profiles for all components (bold line) and the two magnetic compo- 
nents separat e ly (tr iangles and boxes) are presented together with models of 
iMaltbv et al.l (|1986[ ) (thin solid lines, models E and L from top to bottom). 
An indication of the depths probed by our method is given below where the 
normalized line depression contribution functions (CF) are plotted for Fe I 
(dash triple dotted), TiO (dashed), and MgH (dot dashed) lines. 



5 Conclusions 

The aim of this study is to improve a sunspot model atmosphere. Simultaneous 
full Stokes inversions were carried out employing atomic and molecular lines such 
as Fe I, TiO, and MgH. This allowed us to increase the range of temperature 
sensitivity and improve the accuracy in depths ranging from lo g; r = —4 to 
logr <0. As compared to existing umbral temperature profiles (IMaltbv et al.l 
1986I ), the sunspot we analyzed can be best modeled with a warmer atmosphere 
consisting of three components accounting for unresolved umbral structures. 

A detailed disc ussion of this sunspot model atmosphere is in preparation 
f|Wenzel et alJl2010h . 
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